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i-^' The B-meson light-cone distribution ampHtude (LCDA) is defined as the matrix element of a 

quark-antiquark bilocal light-cone operator in the heavy-quark effective theory (HQET) and is 
, a building block of QCD factorization formula for exclusive B-meson decays. When the corre- 

0^ ! sponding bilocal HQET operator has a light-like distance t between the quark and antiquark fields, 

I the scale ^ l/t separates the UV and IR regions, which induce the cusp singularity in radiative 

CsJ ' corrections and the mixing of multiparticle states in nonperturbative corrections, respectively. We 

lO ' treat the bilocal HQET operator based on the operator product expansion (OPE), disentangling 

' the singularities from the IR and UV regions systematically. The matching at the next-to-leading 

■ order is performed in the MS scheme with a complete set of local operators of dimension 

d <5, through a manifestly gauge-invariant calculation organizing all contributions in the coordi- 
, nate space. The result exhibits the Wilson coefficients with Sudakov-type double logarithms and 

H ' the higher-dimensional operators with additional gluons. This OPE yields the B-meson LCDA 

for t less than ^ 1 GeV , in terms of A = — and the two additional HQET parameters as 
matrix elements of dimension-5 operators. The impact of these novel HQET parameters on the 
integral relevant to exclusive B decays, Xb, is also discussed. 
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For the exclusive S-meson decays, such as B ^ nn, P7, . . ., systematic methods have been de- 
veloped using QCD factorization based on the heavy-quark limit In the corresponding fac- 
torization formula of the decay amplitude, essential roles are played by the light-cone distribution 
amplitudes (LCDAs) for the participating mesons, which include nonperturbative long-distance 
contributions. In particular, in addition to the LCDAs for the light mesons n,p, etc., produced in 
the final state, the LCDA 0+ for the B meson, defined as the vacuum-to-meson matrix element [Q], 

hit,^) = ^(0|^-(m)Pe'^^o^A„.A(A.)^^^;^^(0)|^(^)^ ^ J dcoe-""'Mco,^) , (1) 

also participates in processes where large momentum is transferred to the soft spectator quark via 
hard gluon exchange [|l|-|3|]. Here, the bilocal operator is built of the ft-quark and light-antiquark 
fields, /iv(0) and q{tn), linked by the Wilson line at a light-like separation tn, with as the light- 
like vector (n^ = 0, n-v = 1), and representing the 4-velocity of the B meson; a difference 
between (|l|) and the familiar pion-LCDA is that hy (0) is an effective field in the heavy-quark effec- 
tive theory (HQET). jj. denotes the scale where the operator is renormalized, and is the decay 
constant in HQET, = —i{0\qAY5hy\B{v)). The RHS in (|l]) defines the momentum representa- 
tion, with G)v+ denoting the LC component of the momentum of the light antiquark. 

The "IR structure" of ([l|), studied using constraints from the equations of motion (EOM) and 
heavy-quark symmetry [^, as well as the "UV structure", calculated in the 1-loop renormalization 
of the bilocal operator in (|l|) [^, is notoriously peculiar compared with the pion LCDA. For a 
full description of (|I|) which would involve a complicated mixture of the IR and UV structures, 
we first calculate the radiative corrections, taking into account hard and soft/coUinear loops. The 
one-particle-irreducible 1-loop diagrams (ILDs) for the 2-point function {q{tn)AYshy{0)) of (|l]) 
yield 0] ((•••) = (01 ••• |B(v)), the Wilson hne is suppressed, and Cp = {N^ - l)/i2Nc)) 



lLDs=^ f'd^ 
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in D = 4 — 2e dimensions and Feynman gauge, where L = In [i{t — /0)/xe^^] with the MS scale /i 
and the Euler constant Ye- The "vertex-type" correction that connects the light-like Wilson line 
and q{tn) in ([T]) is associated with only the massless degrees of freedom and yields the scaleless 
loop-integral that gives the term with the "canceling" UV and IR poles, 1 / Suv — ^/£ir, and with the 
"plus"-distribution (c^ /(I — ))+ as the splitting function; this term is identical to the corresponding 
correction for the case of the pion LCDA. The other terms in (§) have "non-canceling" UV and IR 
poles: another vertex-type correction around a "cusp" between the two Wilson lines, the light-like 
Wilson line of (|l|) and the time-like Wilson line from /Jv(0) = Pexp[/g J^^dXv ■ A{Xv)]hy{—oov), 
gives the terms proportional to 5(1 — i^), which contain the double as well as single UV pole, 
corresponding to the cusp singularity [^. The "ladder-type" correction, connecting the two quark 
fields in @, gives all the remaining terms in (^), which contain the IR poles and are associated 
with not only the bilocal operator in (^, but also the higher dimensional operators; the ellipses in 
(^) are expressed by the operators involving two or more additional covariant derivatives. 

The renormalized LCDA is obtained by subtracting the UV poles from (Q) with the trivial 
quark self-energy corrections complemented. Here, the term with the plus-distribution (i§/(l — 
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1^))+ is analytic (Taylor expandable) att = 0, similar to the pion LCDA, but the other terms are not 
analytic due to the presence of logarithms L, [§, |8]]. In particular, the nontrivial dependence of 
the latter terms on f/i through L implies that the scale l/t separates the UV and IR regions. Thus, 
we have to use the operator product expansion (OPE) to treat the different UV and IR behaviors 
simultaneously: the coefficient functions absorb all the singular logarithms, while, for the local 
operators to absorb the IR poles, we have to take into account many higher dimensional operators. 
Such OPE with local operators is useful when the separation t is less than the typical distance scale 
of quantum fluctuation, i.e., when t < I//1. We note that an OPE for the B-meson LCDA (|I|) was 
discussed in [^], taking into account the local operators of dimension d <4 and the NLO {0{as)) 
corrections to the corresponding Wilson coefficients in a "cutoff scheme", where an additional 
momentum cutoff Aj/y (» Aqcd) was introduced, and the OPE, in powers of 1 /Ayy , was derived 
for the regularized moments, Mj = f^^"^ d(0(oj(j)+{(0,lJ,), in particular, for the first two moments 
with j = 0, 1; note, Mj ^ 00 as Aj/y ^ 0° [§]. Here, we derive the OPE for (|l|), taking into account 
the local operators of dimension d < 5 and calculating the corresponding Wilson coefficients at 
NLO accuracy. Following the discussion above, we carry out the calculation for f < l/^u in the 
coordinate space and in the MS scheme, so that there is no need to introduce any additional cutoff. 

The most complicated task is the reorganization of contributions from (many) Feynman dia- 
grams in terms of the matrix element of gauge-invariant operators including higher dimensional 
operators, in particular, the three-body operators of dimension 5, such as qGapikj^hy, with the field 
strength tensor Gap [@, §]. To derive the NLO Wilson coefficients associated with such operators, 
we have to compute the 1-loop diagrams for the 3-point function, as well as those for the 2-point 
function as in (^, where the former diagrams are obtained by attaching the external gluon line to 
the latter diagrams in all possible ways. We employ the background field method [[Io|], where the 
background fields represent the nonperturbative long-distance degrees of freedom and satisfy the 

(c) 

exact classical EOM. We use the Fock-Schwinger gauge, x^A^i (x) = 0, for the background gluon 
field A^^^ . This gauge condition is solved to give a[[' (x) = Jq duux^ G^^j^ (ux) | [T0[ ] , which allows us 
to reexpress each Feynman diagram in terms of the matrix element of the operators associated with 
the field strength tensor. Also, this ensures that the Wilson line in (|l|), as well as the heavy-quark 
propagator, does not couple directly to the background gluons while a massless quark or gluon 
propagator couples to them. With the matching in the MS scheme, we obtain M\ the OPE, 



k=l k=\ 

(3) 

where the summation is over a basis of local operators of dimension-J, (k = 1,2, . . .), defined 

as ^ = q^Ysh, {^f ^} = {q{in ■ 'D)4ysh,, q{iv ■ 'B)iiysK}, and {fff} = {q{in ■ ^DfijiK, 
q{iv'D){in-'D)4y5h,., qiiv-'Dfiyshy, qigGapv^'nPj^Ysh^, qigGapfriP^jsh,., qigGapfv^r^JiK, 
^gGap'^'^^^Ys^v}, with another light-like vector, ff- = 0, as = {n^ + n^)/2. The NLO Wilson 
coefficients are obtained as 

Cf\t,,) =1-^(2L^ + 2L + ^) , Ci^)(.,,) = -.[1-^(2L^ + L + ^) 
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(4) 



3 



OPE for B-meson distribution amplitude and dimension-5 HQET operators 



Kazuhiro Tanaka 



and, for the explicit form of C^\t,lJ.),cl^\t,iJ.),. . . ,cf\t,}JL), we refer the readers to Here 
and below, /i is the MS scale, and = av(/i). The double logarithm in the coefficient functions 
originates from the cusp singularity (see (^). The 1-loop corrections for the 2-point function induce 
all of the above ten operators using the EOM, while those for the 3-point function induce only 
'^4^5 6 7 associated with the field-strength tensor; as a result, the coefficients cfl^-i(t^}x) involve 
the terms proportional to the color factor Co = Nc as well as to Cf 

Taking the matrix element (• • •) = (0| ■ • • |B(v)) of (|3|), we can derive the OPE form of the 
B-meson LCDA (|I]). The matrix elements of the local operators in (|3|) are known to be related 
to a few nonperturbative parameters in the HQET, using the EOM and heavy-quark symmetry as 
demonstrated in [|, |]]: (^f ') = AiF{pL)K/3^, {ff!^^) = iF{ix)K, with F of % and A = mB-ntb, 
representing the mass difference between the B-meson and Zj-quark, and all seven matrix elements 
for the dimension-5 operators can be expressed by F, A and two additional HQET param- 
eters Xe and A//, which are associated with the chromoelectric and chromomagnetic fields inside 
the B meson as {qgE ■ aysh^,) = and {qgH ■ aysK) = iF {ix)Xjj{}x) , respectively, in the 

rest frame where v = (1,0). As a result, we obtain |0] the OPE form for the LCDA (I), 
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which takes into account the Wilson coefficients to 0{as) and a complete set of the local operators 
of dimension d < 5. Fourier transforming to the momentum representation and taking the first two 
(j = 0, 1) regularized-moments, Mj = f^"^ d(0(0^ (^+{(0,pi), the contributions from the first line in 
(§), associated with matrix elements of the dimension-3 and -4 operators, coincide completely with 
the result obtained in The second and third lines in (^ are generated from the dimension-5 
operators. Our OPE result (^ "merges" the UV and IR structures peculiar to the B-meson 
LCDA, so that it embodies novel behaviors that are completely different from those of the pion 
LCDA: /I and t are strongly correlated due to the logarithmic contributions, L = In [i{t — /0)/xe^'^], 
from radiative corrections, so that the DA is not Taylor expandable about f = 0, which in turn 



implies the UV divergence in the moments 



as A, 



UV 



oo. The DA receives 



the contributions from (many) higher dimensional operators, in particular, from those associated 
with the long-distance gluon fields inside the B-meson. It is instructive to draw a comparison 
with the previous results, concerning UV or IR structure: one can prove [Q] that (§) satisfies the 
renormaUzation group equation for (|l]), which is governed by the evolution kernel determined by 
the (single) UV poles in (^. On the other hand, (^ reveals that the solution of the EOM constraints 
for ([l|), which was obtained in ||5|], is subject to additional effects from radiative corrections, see 
[^] for the detail (see also [||]). Such corrections to the EOM constraints at order in perturbation 
theory is peculiar to the heavy-meson LCDAs in the HQET and does not arise for the case of the 
(higher twist) LCDAs for the light mesons, 71, p, etc. [11]. 

Our OPE form (|5|) allows us to parameterize all nonperturbative contributions in the B-meson 
LCDA ([I]) by three HQET parameters. A, Xe and A//, and gives a model-independent description 
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• • • • i;c = . 6 GeV" 

t:c = 1.0 GeV" 

Zc=1.4 GeV- 

-■■ Eq.(6; 




Figure 1: The B-meson LCD A at /i = 1 GeV using the OPE (left) and its continuation with a model (right). 



of the B-meson LCDA when f<l//x(<l /Aqcd), taking into account the UV and IR structures si- 
multaneously. Here, we evaluate (||) at the scale jJ. = 1 GeV: A = — m/, in (^ is defined by the b- 
quark pole mass m/,. Following [Q, we eliminate A in favor of a short-distance parameter, Ada, free 
from IRrenormalon ambiguities and written as A = A£)^(/i) [1 + (7/167r)Cf a^] — (9/87r)/xC/r a,, to 
one-loop accuracy; Aoa(m) can be related to another short-distance mass parameter whose value 
is extracted from analysis of the spectra in inclusive decays B —>■ X^y and B —>■ XJ V, leading to 
Az)a(m = 1 GeV) ~ 0.52 GeV [^]. For the novel parameters associated with the dimension-5 op- 
erators, we use the central values of = 0.11 ±0.06 GeV^ = 0.18 ±0.07 GeV^, at 
/I = 1 GeV, which were obtained by QCD sum rules [Q]; no other estimate exists for or A//. We 
calculate ^ for imaginary LC separation, performing the Wick rotation t —^ —ix [Q, ^. 

The results for 0+(— /t,/x = 1 GeV) using (||) are shown as a function of T in the LHS of 
Fig. 1 [^]: the wide-soUd curve shows the whole contributions of (|5|), while the narrow-solid curve 
shows the result for a, — > 0; the NLO perturbative corrections are at the 10-30% level for moderate 
T of order 1 GeV^' ~ 1//I, while they are very large for T ^ because of singular logarithms 
and L. The dashed and dot-dashed curves show the contributions of the first two terms and the 
first line in (||), respectively, associated with the operators of dimension d = 3 and d <4, while the 
dotted curve gives the results of (|5|) when Xe = Xh = 0. For moderate T, the contributions from 
the dimension-4 operators suppress the DA by 30-40%, but the dimension-5 operators, in contrast, 
lead to enhancement by 10-20% with significant effects from Xe and Xh- Our B-meson LCDA (^ 
indeed works up to moderate LC distances T, where the hierarchy among the dashed, dot-dashed, 
and wide-solid curves demonstrates convergence of the OPE (^. 

The two-dot-dashed curve in the LHS of Fig. 1 shows the behavior of the two-component 
ansatz by Lee and Neubert [P], which is given in momentum space as 



0t^(«,At) 



6{(0-(0t) 



71(0 



2 n 



+ 



4Am 
3(0 



In^ 



(6) 



where the second term reproduces the correct asymptotic behavior of the DA ([T]) for ft) » Aqcd and 
the first term represents the nonperturbative component modeled by an exponential form [Q], with 
cOt = 2.33 GeV, = 0.963, and ft^o = 0.438 GeV at = 1 GeV; these parameters are fixed by match- 
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A| - 


= 0.11 GeV2 


= 0.18 GeV2 




AJ = 


A2=0 


[GeV-1] 




(Oo [GeV] 


' [GeV-i] 




(Oo [GeV] 


1 [GeV-i] 


0.4 


0.816 


0.257 


3.11 (0.23 + 2.88) 


0.832 


0.301 


2.69 (0.23 + 2.46) 


0.6 


0.850 


0.306 


2.70 (0.35 + 2.35) 


0.899 


0.394 


2.19 (0.35 + 1.84) 


0.8 


0.852 


0.308 


2.69 (0.47 + 2.22) 


0.966 


0.461 


1.99 (0.46+1.53) 


1.0 


0.858 


0.313 


2.66 (0.58 + 2.08) 


1.11 


0.572 


1.79 (0.56 + 1.23) 


1.2 


0.910 


0.349 


2.51 (0.67 + 1.84) 


1.55 


0.839 


1.56 (0.64 + 0.92) 


1.4 


1.09 


0.456 


2.22 (0.76+1.46) 


4.43 


1.95 


1.32 (0.71+0.61) 


1.6 


1.81 


0.777 


1.87 (0.83 + 1.04) 


9.82 


-4.55 


1.11 (0.77 + 0.34) 



Table 1: Parameters of the model function N/ {tcOq + 1)^ for T > Tc with different values of T^, and the 
results of the inverse moment A^'(/x) at /i = 1 GeV, with the first and second numbers in the parentheses 
denoting the contributions from the first and the second terms in the RHS of (7). 



ing the first two (j = 0, 1) cut-moments J^"^ d(0(0^ ^^{(0,jx) with the OPE for the corresponding 
cut-moments Mo,i derived in where the operators of dimension d < 4 and the corresponding 
Wilson coefficients at NLO are taken into account. For T < 1 GeV^^ the Lee-Neubert ansatz 
shows behavior similar to with A^ = A// = substituted; note that the first term of (Q) produces 
particular contributions associated with the operators of dimension-5 and higher. 

For T S> 1 GeV^\ the contributions associated with higher-dimensional operators become 
important, and the OPE diverges (see (||) and Fig. 1); thus, one has to rely on a certain model 
for the large T behavior and comiect the model-independent descriptions at small and moder- 
ate T to that model. The results in Fig. 1 suggest the possibility of connecting the behavior for 
T < Tc (Tc ^ 1 GeV^^) given by our OPE form (||) to that for t > T^, given by the coordinate- 
space representation of the first term of (||), J^dcoe^'"'^ [Nco/cOq) e^'^l'^ = N / {z(Oq + 1)^. Here, 
A'^ and (Oq can be determined such that both the resulting total DA 0+(— /t,/x) and its derivative 
50+(— /t,/i)/5t are continuous at T = Tc. In the LHS of Table 1, we show the values of A'^ 
and (Oq obtained by solving the corresponding conditions of the continuity for = 1 GeV. (The 
RHS of Table 1 shows the results that would be obtained by solving the similar continuity con- 
ditions with As = A// = 0.) In the RHS of Fig. 1, the wide-solid and two-dot-dashed curves are 
same as those in the LHS, and the dotted, solid-gray, and dashed curves show the behavior of the 
above model function N / {t(Oq + 1)^ for T > Tc with = 0.6, 1.0, and 1.4 GeV^\ respectively, 
using the corresponding values of A'^ and (Oq in the LHS of Table 1 ; these three curves behave as 
^ N/{(O^T^) at large T, with larger A'^/Wq than those of @ and the RHS of Table 1. Indeed, we can 
show that N/(0^ = (9/4A2^) { 1 + t.Ada [^i/^)A + ^«/(2A|,a) - l] } + • • using the continuity 
of 0+(— /t,/i), d^^{—iT,lx)/dT at T = Tc, and thus the contributions of Af and Xh enhance N/(Oq. 

Using these results, we calculate the first inverse moment of the LCDA, 

Xs\^)=l d(o ^+^^'^^ = r dx^+{-iT,^)+ I dz^+i-iz,^), (7) 

JO (0 Jo Jz, 

which is of particular interest for the QCD description of exclusive S-meson decays. We substi- 
tute and the model function, N / {x(0o + 1)^, into the first and the second terms in the RHS, 
respectively, and the results are shown in Table 1 for each value of [^. The "stable" behav- 
ior observed for 0.6 GeV"^ < < 1 GeV"' in the LHS of Table 1 and in the RHS of Fig. 1 
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suggests that X^^n = 1 GeV) ~ 2.7 GeV^^ i.e., Xb{ix = 1 GeV) ~ 0.37 GeV. This value of Xb 
is somewhat smaller than the previous estimates that include nonperturbative and/or perturbative 



QCD corrections [|, |, (e.g., (|) gives = 1 GeV) ~ 0.48 GeV). A value of Xb that is as 



small as our value was adopted in [|l|]. Note that in the RHS of Table 1 with Xe.h = 0, the stable 
behavior is not seen as clearly as in the LHS, and assumes larger values than in the latter. These 
results demonstrate that the novel HQET parameters, Xe and Xh, associated with the dimension-5 
quark-antiquark-gluon operators, could lead to smaller value of Xb- In particular, using the values 
A| = 0.17 GeV^ X^ = 0.25 GeV^, which correspond to their upper bound from the QCD sum rule 
estimate at = 1 GeV [Q], we find that the wide-solid curve in Fig. 1 becomes further enhanced in 
the moderate T region, so that gives Xb{ij. = 1 GeV) 0.2 GeV or smaller. 

To summarize, we have derived the OPE that embodies both the notorious UV and IR be- 
haviors of the B-meson LCDA, including all contributions from the local operators of dimension 
d <5 and the corresponding Wilson coefficients at NLO accuracy. This OPE provides us with the 
most accurate description of the B-meson LCDA for distances less than ~ 1 /Aqcd- We have also 
used the model-independent behaviors from our OPE to constrain the long-distance behavior of 
the LCDA and estimate the first inverse moment X^^ . The results demonstrated the impact of the 
novel HQET parameters, associated with the dimension-5 quark-antiquark-gluon operators. 
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